Ratcheting Heat Flux against a Thermal Bias 
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Merely rocking the temperature in one heat bath can direct a steady heat flux from cold to 
hot against a (time-averaged) non-zero thermal bias in stylized nonlinear lattice junctions that are 
sandwiched between two heat baths. Likewise, for an average zero-temperature difference between 
the two contacts a net, ratchet-like heat flux emerges. Computer simulations show that this very 
heat flux can be manipulated and even reversed by suitably tailoring the frequency (< 100 MHz) of 
the alternating temperature field. 
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The generation of heat flow and its controlled manip- 
ulation presents an ever-growing endeavor for mankind. 
In quest of its technological solution we could witness 
substantial progress over the least decades, with first se- 
rious efforts beiiig achieved that can be traced back to 
the early 1960's 0, S, i, S, H ■ This underlying challenge 
does not present plain sailing because phonons arc by 
far more difficult to control than electrons and photons. 
The recent years, however, have given headway to new 
advances. In particular, thermal rectifiers have been de- 
signed theoretically [1, 0, S, M, [13, 11 1 with a first exper- 
imental realization put forward with help of asymmetric 
nanotubes Moreover, using the negative differential 
thermal resistance , a thermal transistor has been pro- 
posed [13] , which is able to control heat flow much like a 
Field-Effect- Transistor(FET) docs for electric currents. 
Even different thermal logic gates [l3| have been con- 
ceived. All this progress implies that phonons, - tradi- 
tionally being regarded rather as a nuisance -, can in fact 
be put to work constructively in order to carry and pro- 
cess information effectively. Altogether, this has giving 
cause for a new discipline - phononics -, i.e. the science 
and engineering of phonons [15| . 

In addressing this theme let us remind again of the 
original formulation of the second Law by Rudolf Clau- 
sius in 1850 which states that heat cannot spontaneously 
flow from a subsystem at lower temperature to a cou- 
pled subsystem at higher temperature^. Thus, in order 
to generate a steady heat flow against a thermal bias, 
or even generate heat flow in absence of a thermal bias, 
the system necessarily must operate away from thermal 
equilibrium, beyond the limiting realm of the 2nd Law. 
A typical such situation that comes to mind is the Peltier 
effect where a steady heat flow is generated via imposing 



a stationary electric current across an isothermal junc- 
tion of two different materials. 



With this work we propose via computer simulations 
an intriguing new scheme that does not require the re- 
source of a stationary non-equilibrium bias in the form 
of e.g. a stationary electric current but instead com- 
bines the elements of an asymmetric lattice structure 
with a non-biased, temporally alternating bath tempera- 
ture. Dwelling on ideas from the field of Brownian mo- 
tors [TtI . [isl . [igj . - originally devised for Brownian 
particle transport, - we here attempt to direct a priori 
energy (heat) across a spatially extended nonlinear lat- 
tice, see Fig. [1] against an external thermal bias. This 
objective is thus similar in spirit for devising machines 
and devices that can pump heat on a molecular scale 
[13, liSl • doing so, the lattice system is brought 

into contact with two heat baths, with one bath sub- 
jected to a time-varying temperature. Taken alone, this 
nonlinear lattice structure exhibits a thermal diode ef- 
fect 0, 0] which can be exploited to function as a heat 
ratchet device when an additional source of nonequilib- 
rium, - here realized with a rocking bath temperature - , 
is present. Then, merely rocking the temperature in one 
heat bath can induce dynamically a finite bias between 
the two heat baths, being held at the same time-averaged 
temperature. This novel nonequilibrium ratchet feature 
can be utilized (i) to reverse the flux, to (ii) direct heat 
flow from cold to hot against an average thermal bias, 
or even to (iii) turn a regime with a negative differential 
thermal resistance (NDTR) into a regime with positive 
DTR, and vice versa. 



[1] The correct formulation of the second law involves quantities at 
(constrained) thermal equilibrium only; in particular, no time 
variable t enters the formulation of the 2nd Law |26| 



Explicitly, we study numerically a system composed of 
a Frenkel-Kontorova (FK) nonlinear lattice of period a 
which is weakly coupled to a Harmonic lattice (HL), each 
consisting of N/2 atoms of identical mass m. This setup 
is shown in Fig. [H with the FK lattice on the left and 
Harmonic lattice on the right. The FK-HL conflguration 
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FIG. 1: (Color-online) Schematic setup of a weakly coupled 
Frenkel-Kontorova lattice to a Harmonic lattice, being cou- 
pled to two heat baths at temperatures Tj^it) := Tl and T^. 



is governed by the Hamiltonian: 
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Herein, qi = Xi — ia denotes the displacement from 
equilibrium position ia for i-th atom, a is the lattice pe- 
riod, fci and Vl are the spring constant and the on-site 
potential of the FK lattice, kint is the coupling strength 
between the FK and the Harmonic lattice, and kn is the 
spring constant of the Harmonic lattice. Fixed bound- 
ary conditions, yielding ~ Qn+i = 0, have been 
employed. The 1-st atom and the iV-th atom are put 
into contact with two Langevin heat baths possessing 
temperature T^it) := Tl and T^, respectively. Gaus- 
sian white noise are used, namely, (^i/Ar(i)) = and 
(^i/Af(i)Ci/Ar(0)} = 2kBr]Ti^/Rd{t). ks is the Boltzmann 
constant and rj denotes the coupling strength between 
system and heat bath. The time-varying heat bath tem- 
perature T^lt), oscillates dichotomously at angular fre- 
quency Lo and driving strength A. The used bath tem- 
peratures thus read explicitly: 

mt) := Tl = To{l + A + A- sgn{sin ujt)) , 

Tr = To(l - A) , (2) 



where Tq = {TL{t) + Tfl)/2 is the temporally av- 
eraged environmental reference temperature, 2A = 
{TLit) —Tjij/Tf) denotes the normalized temperature dif- 
ference, and sgn(sina;t) provides the dichotomous, time- 
dependent temperature variation. The time scale uo of the 
temperature manipulation of the heat bath is assumed 
to vary much slower than the time scale Tieq to reach 
local thermal equilibrium; i.e. a; << t^"^. This time 
scale for good thermal conductors is typically a function 
of temperature; it is of the order of the time scale of 
the electron-phonon relaxation time which normally de- 
creases with decreasing temperature. For good metals 
this time scale is of the order of 0.1 — 1 picoseconds. 

We next use dimensionless parameters by measuring 
positions in units of a, momenta in units of [a(mfci)^/^]. 
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FIG. 2: (Color Online) The numerically (10® realizations) 
evaluated, periodically varying asymptotic local heat flux 
Ji(t) over the full driving period 2tt/uj — 100 is compared 
with the rocking bath temperature TL{t), at driving strength 
A = 0.5 and zero bias A = for To = 0.09. The lattice length 
is A'' = 50 -I- 50. The numerically cycle averaged, spatially 
homogenous heat flux equals J = —4.92 ■ 10^® (i = 5) and 
J : 



-4.93 ■ 10"" (j = 10). 



temperature in units of [kLO?' /ks], spring constants in 
units of k^, frequencies in units of [(m/fci)^/^] and en- 
ergies in units of [fc^a^]. In particular, we set in our 
simulations fc^ = 1,1/^ = 5,fci„t = 0.05, fc/j = 0.2. For a 
typical situation, the dimensionless temperature is set at 
To = 0.09. This yields with kLa^/ke - 10^ A' - lO^K a 
physical temperature of the order Tq ~ 90K — 900K. The 
equations of motion are integrated by the symplectic ve- 
locity Verlet algorithm with a small time step h = 0.005. 
The system is simulated for a total time tt = 2 - 10^. The 
chosen optimal coupling strength of the heat bath is fixed 
at ?7 ~ 0.5. 

The asymptotic heat flux Ji{t) is assuming the pe- 
riodicity of the external driving period 2tt/lu after the 
transients have died out. This fact is assured for all 
of our chosen frequencies uj after a simulation time of 
st ~ O(IO^). At those asymptotic long times t the heat 
flux equals the noise average Ji{t) — ki {qi{qi — qi-i)) 
where ki = k^ for i = 2, . . . , N/2 and ki = ku for 
i = {N/2) + 2, . . . , N , being here evaluated in the com- 
monly employed way, cf. in Refs. U, 13, 24 1. 
The static heat flux J then follows as the cycle average 
over a full temporal period: J = Jo'^^'^ Ji{t)dt, which 
with ergodicity being valid equals as well the long time 
average, i.e. J = kLq2{q2 - qi) = kaqiyiqN - qN-i)- In 
Fig. [2] we depict the resulting periodic variation of the 
heat flux Ji{t) vs. the externally applied temperature 
variation TL{t). This nonlinear response exhibits a char- 
acteristic phase lag vs. the perturbation T^it) and is 
dynamically biased to yield a nonzero temporal average. 

With the introduction of a temporally alternating tem- 
perature field in Eq. ([2]), we thus achieve a controllable 
manipulation of heat flow by an externally adjustable pa- 
rameter, i.e. the driving frequency lo. The static thermal 
bias A has been set to zero. The only resource driv- 
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FIG. 3: (Color-online) Time averaged heat flux J vs. the 
angular driving frequency uj for different lattice lengths A*' = 
50 + 50 and AT = 100 + 100 and To = 0.09. The two arrows 
mark the heat flux calculated in the adiabatic limit. 



ing heat across the junction thus is the non-equilibrium, 
alternating temperature field Tl (t) , which generates pos- 
itive and negative temperature variations in the first and 
second half of driving period. As a result, the direc- 
tion of heat flow will tend to reverse each half driving 
period. In the adiabatic limit; i.e. a; 0, the alter- 
nating temperature Tl (t) can be expressed by two oppo- 
site static thermal bias values, yielding the average heat 
current approximately as the averaged heat current for 
two opposite static thermal bias values, see the two hor- 
izontal arrows depicted in Fig. [31 In contrast, in the 
fast-driving limit cu oo, the left-end atom will expe- 
rience a time-averaged constant temperature. This cor- 
responds to thermal equilibrium, yielding J when 
UJ oo. In Fig. [31 the numerically determined aver- 
age heat current is depicted as a function of the driving 
frequency to for the driving amplitude A — 0.5. In full 
agreement with our predictions, a finite heat current J 
emerges in the adiabatic limit w ^ 0, becomes dimin- 
ished in the non-adiabatic limit and essentially vanishes 
for large cu. At adiabatic driving the values of J agree 
well with the numerical values determined from the adi- 
abatic approximation. A tantalizing observation is that 
J does not vanish monotonically. Remarkably, at some 
intermediate value lo, the heat flow J crosses zero and 
subsequently reverses its direction upon further increas- 
ing UJ. Consequently, the direction of net heat flow can 
be manipulated by suitably tailoring the frequency of the 
temperature variations. 

This interesting reversal of the heat flux can be re- 
lated to the thermal response time of the system. The 
non-rocked FK lattice obeys Fourier's law [2^. Thus, 
the corresponding temperature variations T(t, x) obey 
the diffusion equation: dT/dt = {n / c)d'^T / dx^ ^ where n 
denotes the heat conductivity and c is the specific heat. 
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FIG. 4: (Color-online) (a) Stationary heat flux J vs. driving 
frequency lo for different lattice lengths A'^ = 50 -f 50 and 
iV = 100 + 100 at To = 0.09 and a driving strength A = 0.5. 
(b) The effective temperature proflles of four selected points 
in (a) at length iV = 100 -I- 100. 



The solution follows a Gaussian wave packet T{t, x) — 
l/2^^T:Kt/cexjp{~x'^c/4,Kt). The thermal response time 
can now be estimated as the time span for the energy to 
diffuse across the system, i.e. r ~ cN"^ /An. At tempera- 
ture To ~ 0.09, the FK lattice assumes the numerical val- 
ues K ~ 0.5 and c ~ 1. Thus, the characteristic frequency 
scale ujc of the system can be estimated as Wc = 27r/T. 
This characteristic frequency scales inversely with the 
system size oc iV~^. For N = 100 we then find lOc ~ 10""^ 
and a roughly four times smaller value for N = 200. 
These two estimates for ujc are in good agreement with 
those observed numerically in Fig. [31 Taking the physical 
unit of frequency, i.e. wq = [(m/fc^)-^/^] ~ lO^'^sec^^ into 
account, lOc ^ 10~^ corresponds to a typical physical fre- 
quency for microwaves of ~ 10^° sec~^. The theoretically 
predicted red shift oc 7V~^ for ujc with increasing system 
size is nicely corroborated by our numerical results. 

To gain additional insight into this reversal of heat flow 
we investigate the local temperature variations across 



4 



o. \ 
-o \ 

"q \ 


FK+FK 

N = 50+50, A = 
— • — A = 0.5 

--o-A = 0.4 
9 A = 0.2 


bk 

»-■<»- 'A 

"a.. i 







10^ 



10" 



10^ 



10" 



10" 



10" 



FIG. 5: (Color-online) Stationary heat flux J vs. driving fre- 
quency u! for different driving strengths A of the temperature 
modulation for a size A'' = 50 -I- 50 at To = 0.09 and zero 
thermal static bias A. 



the junction at different driving frequencies to. After 
evolving the system over a long total simulation time 
tt = 2 ■ 10* this local temperature of i-th atom is evalu- 
ated from its temporal Jxjng time average of the kinetic 
energy, i.e. Teff{i) — qf- In doing so, we switch to a 
FK-FK configuration because the harmonic lattice know- 
ingly is not able to build up a temperature gradient 24 1 . 
The employed right sided FK lattice has the parameters 
{kfj = 0.2, Vr = 1). In Fig. [UJa) a similar heat current 
modulation as for the FK-HL configuration is observed 
for the FK-FK junction with A = and A = 0.5. The 
effective temperature profiles of four numerical runs, de- 
noted as (a,b,c,d) in Fig. [D^a), are depicted in Fig. ID^b) 
versus the relative site positions i/N. 

In clear contrast to a non-rocking case (i.e. A = 0) with 
no net temperature bias, a distinct temperature gradient 
now emerges for a rocking temperature T^it). The tem- 
perature profile exhibits a discontinuity at the interface. 
For the situation in (c) where the heat flow reversal J 
approaches its lowest, negative value, the effective tem- 
perature profile becomes rather complicated: Away from 
the bending part (with a negative-valued slope) towards 
the left terminal side the resulting temperature profile 
exhibits an opposite thermal gradient in comparison to 
case (a)-thus indicating that a current reversal occurs. 
This opposite thermal gradient behavior can also be de- 
tected upon comparing the temperature gradients of case 
(b) and (d), exhibiting both a vanishing heat flow J — 0. 
The discontinuity of Teff{i) occurring at the interface 
reaches its maximal value at low frequencies, cf. case (a), 
and increasingly diminishes with increasing driving fre- 
quency beyond the frequency value for reversal, cf. case 
(d). 

We next investigated numerically the role of the driv- 



ing strength A of the temperature modulation versus 
driving frequency u. The results are depicted with Fig. 
dH]). As expected, a lower driving strength consistently 
yields smaller values of the Brownian motor induced heat 
flux J, which vanishes identically when the strength of 
the source of nonequilibrium is approaching zero, i.e. 
A = 0. Because the frequency scale lUc for occurrence 
of heat flux reversal is mainly size dependent oc N~'^ we 
expect a weak dependence of uJc on driving strength A. 
This result is corroborated with Fig. ^ where this char- 
acteristic frequency is practically independent of driving 
strength A. 

Flux-bias characteristics. The flnite ratchet value of 
heat flux J in the absence of static thermal bias A now 
allows for directing heat current against a non-zero ther- 
mal bias A. In Fig. [HI we depict the flux-bias character- 
istics J — A at small driving frequency oj = 1.571 • 10~^ 
(~ 100 MHz), which is in the range of ultrasonic fre- 
quencies. The J — A curve in absence of time-dependent 
manipulation (A = 0) is presented as a reference, obey- 
ing J(A = 0) = 0, in agreement with the 2nd Law. The 
dashed line corresponds to a driving amplitude A = 0.5. 
We now detect a finite heat current J at vanishing ther- 
mal bias A = 0. For negative bias within the range 
A e [—0.23,0], the direction of the heat flow is posi- 
tive. With J taken as the average over a driving period 
this implies that heat flows from cold to hot. The so 
called "stall bias" of the thermal bias that yields J = 
is located around A = —0.23. Moreover, we detect that 
the characteristics of (NDTR) becomes modified as well 
by the switch-on of the alternating bath temperature 
TL{t). The original working range of NDTR at A = 
is A G [—0.6, —0.2]. At A ~ 0.2, this range undergoes a 
shift towards A G [-0.7,-0.3]. For A = 0.5, this very 
NDTR-phenomenon effect can be dynamically eliminated 
all together. This effect is thus of prominent interest for 
the design of a efficient thermal transistor: It allows one 
to change the range of working temperatures of a thermal 
transistor by solely adjusting the strength of the driving 
temperature field. 

In conclusion, we have shown that heat flow across 
a structured nonlinear lattice junction can be efficiently 
controlled by use of a temporally alternating temperature 
fleld. The heat flow becomes directed and even can be 
reversed by suitably selecting the driving frequency. In 
clear contrast to the thermal diode effect, the presented 
Brownian heat rachet physics dynamically generates a fl- 
nite heat flux at zero thermal bias A = 0. Thus, we deal 
with a new phenomenon which, as well, is in distinct 
contrast to the by now common situation of thermally 
assisted, directed particle transport in Brownian motors 
Til 17, 18, 3]. This dynamically induced heat flow can 



be directed against a non-zero, time-averaged net thermal 
bias A. This fact, however, does not necessarily imply an 
active overall cooling of the device, cf. also Fig. 2] (b). 
The size and the shape of NDTR can be manipulated 
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FIG. 6: (Color-online) Heat flux J vs. thermal bias A for 
different driving amplitudes ^ = 0, 0.2, and 0.5. The lattice 
length is = 50-1-50 and To = 0.09. Note that the nonhnear- 
ity in the Frenkel-Kontorva part of the junction is essential to 
obtain the thermal ratchet effect. At large rocking strength 
(A = 0.5) the current-bias characteristics can be manipulated 
to eliminate a NDTR-regime at negative bias values A. 



as well. Nonlinearity, as reflected with the FK-segment, 
is essential for the phenomenon: a junction composed of 
two asymmetric harmonic lattices fails to support a ra- 
chet heat flux. All these phenomena call for beneficial 
applications in the control and management of heat flow 
on the micro- and nano-scale. With a future, more de- 
tailed work [i^l we shall investigate, how the size of the 
directed heat flux can be enlarged by stylizing further 
the nonlinearity of the lattice structure. In particular, 
we expect much larger heat currents when using a Fermi- 
Pasta-Ulam chain containing both, a cubic and a quartic 
interaction potential. 

Our setup constitutes a two-segment rectifier model Q 
and the driving frequency is in a range of (< 100 MHz). 
The finite, directional heat fiux is feasibly detected by 
slowly rocking the bath temperature within the adiabatic 
regime. Given the fact that rectification has been demon- 
strated experimentally in asymmetrically deposited nan- 
otubes we are confident that the appealing heat 

ratchet phenomena presented herein will invigorate ex- 
perimentalists to validate our findings by designing such 
thermal ratchet systems. 
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